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ABSTRACT 

We consider the long term evolution of debris following the tidal disruption of compact stars in 
the context of short gamma ray bursts (SGRBs). The initial encounter impulsively creates a hot, 
dense, neutrino-cooled disk capable of powering the prompt emission. After a long delay, we find 
that powerful winds are launched from the surface of the disk, driven by the recombination of free 
nucleons into a-particles. The associated energy release depletes the mass supply and eventually shuts 
off activity of the central engine. As a result, the luminosity and mass accretion rate deviate from the 
earlier self-similar behavior expected for an isolated ring with efficient cooling. This then enables a 
secondary episode of delayed activity to become prominent as an observable signature, when material 
in the tidal tails produced by the initial encounter returns to the vicinity of the central object. The 
time scale of the new accretion event can reach tens of seconds to minutes, depending on the details 
of the system. The associated energies and time scales are consistent with those occurring in X-ray 
flares. 

Subject headings: accretion, accretion disks — hydrodynamics — gamma rays: bursts 



1. INTRODUCTION 

Gamma Ray Bursts (GRBs) probably involve new- 
born compact objects accreting at rates high enough 
that neutrinos are the dominant source of cooling. Short 
burst s (SGRBs) e xhibit a substantial degree of diver- 
sity (jNakarl l2007f ). although a large fraction of them 
probably come from an old stellar progenitor population 
()Lee fc Ramirez-Rui3l2007t ). while long events (LGRBs) 
have been associated to the collapse of massive stars 
(jWooslev fc Blooml l2006f) . Although two classes are 
probably present, some overlap occurs and making a 
distinctio n for individual ev ents is not entirely straight- 
forward (jZhang et al.l l2009f) . Since the launch of the 
Swift satellite it has been possible to study the onset 
of the afterglow emission, thought to be produced when 
the inertia of the swept up ext ernal matter starts to 
slow down the ejecta appreciablv (IFox fc Meszarosll2006t 
lO'Brien et aH mM iGehrels et all \200Tl ) . ThiThas re- 
vealed in many cases renewed rapid flaring (usually in 
soft X-rays) long after the canonical prompt gamma- 
ray emission has ceased - at timescales that are much 
longer than the usual dynamical or even viscous time 
scales considered for the accretion process itself. In 
SGRBs these flares typically appear 30-60 s after the 
burst and last for tens of seconds as well, containing a 
non-negligible fraction of the radiated energy. There is a 
general consensus that they must involve renewed ac- 
tivity within the central engine, though simple broad 
pulses could arise from an external shock interaction 
jRamirez-Ruiz. Me rloni fc Ree£|[200lH Zhang et al.ll2006l : 
[Nousck et al. 2006). 

We seek to accurately compute the evolution of these 
neutrino cooled structures on long time scales, and con- 
sider mechanisms which might allow for the produc- 
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tion of a new episode of energy release. Disk frag- 
mentation (Perna. Armitag e fc Zhanel [200l l. faU back 
tidal tails (Rosswog 2007 a), strong magnetic fields 



rapidly rotating pulsars (iKluzniak fc RudermanI Il998l : 
iRosswog fc Ramirez-Ruizj 120021: iDai Wang fc Full2006f) 



magnetic h alting (iProga fc Zhanel I2006D an d cooling 
transitions ("Lazza ti. Perna fc BegelmanI |200"8|) , among 
others, may be related to such secondary transients, but 
no definitive answers have been forthcoming. In this Let- 
ter we study the thermodynamical properties of these 
flows in detail in an attempt to determine how they may 
be relevant to their long term evolution, which we com- 
pute numerically. 

2. INPUT PHYSICS AND NUMERICAL STUDY 

Previous work on ^-cooled disks has 
shown the importance of the proper consid- 
eration of thermo dynamic s (IWooslevI 119931: 



MacFadven fc Wooslevir i999: Po pham. \Vooslev fc Frveil 

' 120011: 



19991; LNaravan. Piran fc Kumar 

bi Matteo. Perna fc Nar avan 2002; Kohri fc Mineshi^ 
20021: iBeloborodovl 120031: Setiawan, Ruffert & Janka 
2004 iLee. R amirez-Ru iz fc Page .2005; .Rosswoa ,2007 



Chen fc Be loborodov 20Q7|. The densities and tempera- 



tures are typically p ~ 10 g cm~'^ and T ~ 10 K, and 
we consider an ideal gas of free nucleons and a particles 
in nuclear statistical equilibrium (NSE), black-body 
radiation, relativistic e^ pairs of arbitrary degeneracy, 
and neutrinos, for which we use a simplified two-stream 
app roximation (|L6pez-Camara. Lee fc Ramirez-Ruid 
12009). Neutronization in the optically thick and thin 
regimes is c onsidered by detailed balance of weak 
interactions ()Lee. Ramirez-Ruiz fc Pag^ I2005D . which 
consistently accounts for the optically thick and thin 
regimes. Emissivities are taken from the fitting functions 
given by Ut oh et al. (1996) for pair annihilation and 
iLanganke fc Martinez-Pinedd (|200l[ l for e''^ capture 
onto nucleons, and the energy from photodisintegration 
of a particles is also included in the energy equation. 
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The evolution of the disk is followed with a two di- 
mensional Lagrai iKian Smooth Pa rticle Hydrodynam- 
ics (SPH) code ([Monaghanl |1992| ) in azimuthal sym- 
metry. The central mas s produces a Paczynski-Wiita 
(jPaczyhski fc W iita 1980*) pseudo-Newtonian potential, 
reproducing the position of the last stable orbit for 
a Swarzschild black hole (BH), and accretion is im- 
plemented with an absorbing boundary at rjn = 
2''Sch = 4G'Mbh/c^- We also include a boundary at 
Tout — 200rsch, where outflowing fluid elements are re- 
moved from the calculatio n, and use an a-prescription 
(jShakura fc Sunvaevlll973h for the magnitude of the vis- 
cosity, with 10~^ < a < 10~^. To avoid catastrophic loss 
of resolution a fission routine maintains a predetermined 
minimum number of fluid elements throughout the evo- 
lution. The initial disk mass is 0.003 < Mdisk/A^© < 0.3, 
and the initial BH mass is A/bh — 5Mq. 

3. ISOLATED DISK EVOLUTION AND THE IMPORTANCE 
OF HE SYNTHESIS 

A ring in centrifugal equilibrium at characteristic ra- 
dius R accretes and evolves on a viscous time scale if 
angular momentum is transported outwards and mass 
inwards. The standard approach is to consider that vis- 
cous stresses, and the associated energy dissipation can 
be parameterized through the a prescription, with the 
coefficient of viscosity given by z/ = ac^/flKcp, where Cg 
and f^Kcp are the local sound speed and Keplerian or- 
bital frequency, respectively. The ring spreads radially 
on a timescale ivisc — R/Wv, eventually transferring all 
of the angular momentum to an infinitesimal amount of 
mass at infinity. If the dissipated energy is radiated ef- 
ficiently it remains geometrically thin, with aspect ratio 
h/r ^ Cs/w0 <SC 1. In the context of GRBs, where ac- 
cretion rates can reach ~ 1 Mq s~^, the disk cools by 
neutrino emission and typically the accretion efficiency is 
?7acc = Lu/M(? ~ 0.01 — 0.1. The nature of the mecha- 
nism is quite irrelevant in determining the vertical struc- 
ture of the disk, so long as it is minimally proficien t in 
removing internal energy (|Lee fc Ramirez-Ruidl2002f) . In 
this case the decay in luminosity and accretion rate are 
closely correlated. 

An earlier study considering accretion disks formed im- 
pulsively by the tidal disr uption of main sequen ce stars 
by supermassive BHs by ICannizzo et al.l ()1990l) found 
that the luminosity follows a power law in time, L oc t^^ , 
with (3 ~ 1.2, as does M. A key consideration is that the 
mass of the disk be only remo ved through accretion onto 
the ce ntral object. Recently, iMetzger. Piro fc QuataertI 
(|2008D have carefully computed the evolution of a ring 
of matter in the v-c ooled regime under sim ilar assump- 
tions as the study of IC annizzo et all ()1990 f). The evolu- 
tion displays similar behavior, and they report a decay 
index (3 ~ 4/3. Further, they consider generic solutions 
in which a disk-driven wind is present, extracting mass 
as well as angular momentum from the fluid. In that case 
both L and M enter a phase of rapid decline once mass 
driving through the wind becomes important, typically 
after a fraction of a second. 

Thus motivated, we have computed the long term evo- 
lution of I'-cooled accretion disks under the assumptions 
given in § [21 The early stages [t < 0.5 s) are character- 
ized as shown before by the release of the initial reservoir 
of internal energy left over from the merger phase, and 
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Fig. 1. — Neutrino luminosity for models with Mjisis/M© = 
0.3,0.03,0.003 (top to bottom in each panel) and a = 
0.1,0.01,0.001 (top to bottom panels). The grey lines show the re- 
sult of calculations in which the energy associated with He synthesis 
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are followed by a more gradual decay in which angular 
momentum transport and energy dissipation contribute 
to both draining the disk of matter and powering the 
neutrino luminosity, L^. Depending on the disk's mass, 
the flow can be optically thick to neutrinos in the early 
phases and this can have an effect on the global energy 
release. 

After a longer delay, powerful winds are launched 
from the disk because of He synthesis from the avail- 
able free nucleons, at the characteristic radius where the 
nuclear binding energy of He is roughly equal to the 
gravitational binding energy, GMBH"ip/rwind — 7 MeV. 
When the available mass approaches rwind, it is driven 
away almost isotropically, rapidly depleting the disk and 
shutting down neutrino emission. The effects are seen 
in Figure [1] where Li, is plotted for simulations cover- 
ing two orders of magnitude in initial disk mass and 
strength of viscous transport. For comparison, curves 
of benchmark simulations in which the energy associ- 
ated with He synthesis, 7 MeV/baryon, was not taken 
into account in the energy equation are also shown 
(grey lines), and clearly follow a power law with index 
P 2± 1.3 for a longer tini e , clos e to that predicted by 
iMetzger. Piro fc QuataertI ()2008D (see especially the case 
with a = 0.01). The break to something more akin to ex- 
ponential decay is entirely due to the driving of winds off 
the surface of the disk. The mass flow rates and mechan- 
ical power are (M[10-2 M© s'^], Lwind[10^° erg s'^]) = 
(1.5, 1.3), (0.2,0.15), (0.025,0.01) for initial disk mass 
Mdi.k = 0.3, 0.03, O.OO3M0, respectively, and a ^ 0.01. 

In the reference simulations which do not consider the 
nuclear binding energy, eventually also deviates from 
the simple power law because the mass drains into the 
BH as a result of viscous transport and neutrino cool- 
ing shuts off. As expected, the energy normalization is 
Li, oc Mdisk, since the energy reservoir is gravitational in 
nature, and the temporal breaks scale with the vigor of 
angular momentum transport since a lower a implies a 
longer delay in the mass moving to rwind- 

The thermodynamical conditions clearly affect the evo- 
lution of the accretion disk. Figure shows a snapshot in 
the evolution of one of our models. The two main cooling 
mechanisms are, as mentioned, e^ annihilation and e^ 
capture onto free nucleons. The first requires abundant 
e^ pair creation, while the second only operates if there 
are free nucleons in the fluid. As the transition from a 
particles to a free gas of neutrons and protons is quite 
rapid, and the abundance of pairs decreases exponen- 
tially with degeneracy, these two conditions clearly mark 
the boundaries where the gas is allowed to lie in such a 
diagram in order to cool effectively. If a fluid element 
were to enter a region where cooling is inoperative, the 
associated expansion will quickly cause a drop in den- 
sity and force it back into the cooling region (pressure 
mainly comes from the ideal gas terms in the equation of 
state). As the disk is depleted, it follows the photodis- 
integration line closely until this crosses the degeneracy 
threshold, with T oc p^/"^, and then continues to even 
lower densities along parallel tracks. It is this transition 
to a disk dominated in composition by a particles which 
leads to the production of strong winds and rapidly ex- 
hausts its mass as discussed above. The change occurs 
at logp[g cm-3] ~ 6.5 and T ~ 10^° K. 
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Fig. 2. — Time evolution of a neutrino cooled accretion disk in 
the density-temperature plane. The blue power law marks a line 
of constant degeneracy given by kT oc p^/^. The curved green 
line shows the transition in composition from free nucleons at high 
temperatures and low densities to a particles in the opposite ex- 
treme assuming NSE. The thick black arrow marks the flow of the 
gas as the disk is drained through winds and accretion onto the 
central mass. The spacing between each successive pair of black 
dots along the blue and green curves corresponds to a change in 
the mass cooling rate (in erg g~^ s~^) of one order of magnitude. 
The highest point has the largest value at = 10^'' erg g~^ s~^. 
The red dots show the values of the density and temperature in 
the disk for a snapshot in the evolution of the model with a = 0.01 
and initial disk mass A/^isk = O.OSAfQ. 



4. THE EFFECTS OF MASS INJECTION AT LATE TIMES 

The tidal disruption of a neutron star by a 
more massive BH, or a double neutron star merger, 
ejects stellar material int o long tidal tails through 
the outer Lagrange point (iLattimer fc Schramml [l97l 
iRasio fc Sha"Di73[l9MlLe3l200lHRosswodl2005l . l2007afl ^ 
A fraction of this is has enough energy to escape the 
system, with possibly interesting nucleosynthetic conse- 
quences. The rest is bound to the central object on highly 
eccentric orbits and a range of initial energies. Trajec- 
tories are essentially ballistic and this can provide an 
estimate for the rate of mass fall back. A constant differ- 
ential mass distribution with energy yields homologous 
flow and A/fb = Mo(i/^o)~^^'^, where to and Afo fix the 
normalization, and the particular details depend on the 
equation of state of nuclear matter and the initial mass 
ratio. Simulations indicate that up to 1O~^M0 of ma- 
terial may follow this behavior. The bulk of the mass 
returns to small radii after a few seconds, and the ques- 
tion is whether it can release its gravitational binding 
energy efficiently and at the same time alter the overall 
evolution of the accretion flow. 

We now proceed to describe the evolution of the fall 
back material based on our knowledge of merger event 
properties and the resulting configurations. At some time 
tinjcct — 1 — 10 s, a mass A/fb — 10~^ M© creates a 
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new ring with characteristic radius i?~3 — 5x 10^ cm. 
This is located further out than the original accretion 
torus, typically by a factor of 2-4 because of its some- 
what larger specific angular momentum after ejection in 
the tidal disruption of the neutron star. Note that the 
fall back mass accretion rate is only that at which the 
tidal tails feed the disk, and not the rate at which the 
central object grows. The new ring will dominate the 
evolution of the system provided that the injected mass 
is greater than that remaining in the original disk at 
t = ^inject, Mih{tinjcct) > A^disk (^inject ) ■ For an inte- 
grated value of Mfb ~ 5 x 10~^ Mq we find that this 
is valid at times ranging from 0.4-3 s, depending on the 
initial value of a. The key ingredients are the specific an- 
gular momentum and total mass of the injected material, 
and the strength of viscous transport. Not all combina- 
tions can lead to a clearly defined secondary episode of 
energy release (e.g., if the re-injection is too close to the 
BH and/or too slow, or lacks sufficient mass, the effect 
will be greatly reduced). 

We show the evolution of rings with initial masses 
Mfb = 0.2,0.02 M0 and viscosity a = 0.01,0.001, which 
produce flare-like activity. As the ring spreads radially, 
the temperature and density in the innermost regions rise 
on a time scale that depends on a, and as a result the 
luminosity increases. The fluid is now mostly a parti- 
cles and the main source of cooling is e^ annihilation, 
which is efficient enough to avoid unbinding the gas and 
depleting the disk's mass. The accretion rate and lumi- 
nosity subsequently drop off as a power law with index 
/? ~ 1.3, as in the previous simplified calculations where 
no winds were allowed to form, before making a tran- 
sition to a steeper decay as neutrino cooling becomes 
less efficient (see Figure [3]). The peak neutrino luminos- 
ity reaches log[Ly(erg s"""^)] ~ 48 — 49, and can remain 
above ~ 10'*'' erg s~^ for up to t ~ 100 s. 

Note that we have referred consistently to the release 
of energy through neutrinos in the fiow. This is simply 
because it is the physical ingredient we can accurately 
compute and follow reliably, serving as a tracer of cen- 
tral engine activity. This does not mean that neutrinos 
alone are necessarily capable of powering the observed 
e mission during fiares, as was n oted, e.g. for GRB050724 
bv lFan. Zhang k Progal (|2005( ). Qualitatively, the argu- 
ment applies to any form of energy extraction that is 
dependent upon a mass reservoir, for example, magnet- 
ically driven outflows, with ~ fpc^ in the disk and 
/ < 1. In such a scenario, fields anchored in the disk can 
drive outflows wit h corresponding luminosit ies reaching 
Lb =i lO-^^erg s'^ (|Blandford fc PavndfTOSl) . which are 
high enough to account for the observed power outputs. 

The total neutrino luminosity including the prompt 
phase, plotted in Figure [3l shows that the interval be- 
tween the start of the evolution and the maximum bright- 
ness can reach a minute or more depending on the in- 
tensity of viscous transport. The key ingredient in the 
evolution is that the available mass not be driven off the 
disk before a substantial fraction of it can accrete onto 
the BH and release its gravitational binding energy. The 
computation of the cooling rates is crucial in this respect 
as one needs to maintain an accretion efficiency of a few 
per cent in order to successfully transfer the energy to 
neutrinos. 




Fig. 3. — Neutrino luminosity for the evolution of post merger 
disks and including the contribution from fall back with initial mass 
l^ldisk/ Mq = 0.03 and Mi^/Mq = 0.02. The energy normalization 
scales linearly with the mass. Note that in this plot the time is 
given in seconds. The red and blue curves show the evolution for 
a = 0.01 and a = 0.001, respectively, where the initial decay is 
given by the original disk and the late emission peaking at t ~ 
3, 30 s comes from the material injected at tinjoct — Is. The time 
delay until the flare maximum is clearly dependent on the intensity 
of angular momentum transport. 

The actual shape of the lightcurve is sensitive to how 
the energy is reprocessed into radiation, which in turn 
depends on the Lorentz factor of the outflow before and 
after the mass depletion phase occurs. This would deter- 
mine whether internal dissipation takes place in the sec- 
ondary outflow before it interacts with previously eject ed 
material or not ([Ramirez- Ruiz. Merloni fc ReeslboOll) . 

5. DISCUSSION AND CONCLUSIONS 

The global properties of accretion flows are intimately 
related to the cooling mechanisms present (witness the 
relevance of the two most widely known limiting solu- 
tions, standard Shakura & Sunayev (SS) models and 
Advection Dominated Accretion flows (ADAFs) over the 
past 35 years). For stationary solutions, the mass and 
energy fed into the flow fix the resulting properties. 
However, when temporal variations are allowed, transi- 
tions from one regime to another can naturally occur. 
The solutions exhibited here are well described by the 
self-similar temporal decay c omputed for isolated fluid 
rings with efficient cooling bvlCannizzo et al.] (|199C1|) and 
iMetzger. Piro fc Quataertj (|2008D . but also switch to a 
different mode, and rapidly decay when the disk mass 
drops below the threshold required to sustain adequate 
cooling. In this particular instance, the mechanism re- 
sponsible for mass depletion is a combination of accre- 
tion and neutron and proton recombination into a par- 
ticles. This phase transition then effectively allows the 
secondary episode of accretion to play an important role 
with potentially observable consequences. 

We have shown here how the injection of material from 
tidal tails formed during the initial disruptive encounter 
of compact objects can re-energize the accretion disk pro- 
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vided that the fall-back mass dominates over the remnant 
disk at the time when it is re-injected at its circulariza- 
tion radius. It must be emphasized that the evolution of 
this fluid then procee ds not on the relatively short fall- 
back time scale itself (jRossi &: BegelmanI 12009 ). but on 
the viscous time of the resulting ring. Due to the dy- 
namics of the encounter, the specific angular momentum 
can be large enough to increase the characteristic orbital 
radius and thus account for longer time scales. Upon re- 
turn to the vicinity of the BH, the gas can efficiently cool 
and accrete. The resulting luminosity, integrated energy 
and time scale for delay can account for the observed ac- 
tivity at late times in SGRBs, provided a sufficiently ef- 
ficient mechanism is available to tap the available energy 
(e.g., magneti c fields). The d etails of fall back depend on 
many factors (jLee fc Ramirez -Ruiz 2007), including the 
neutron star equation of state, the initial mass ratio of 
the merging binary, the nature of the components (dou- 
ble NS binary vs. NS-BH system) and the dynamical 
interaction (merger in a binary vs. collision in a dense 
stellar environment). All of these would contribute to 
substantial diversity in the outcome, and it is possible 
that many events would have no observable signal aris- 
ing from such tails, just as not all mergers may produce 
a SGRB. 

Late flares have been observed to date in both long 



Beloborodov, A. M. 2003, ApJ, 588, 931 

Blandford, R.D., Payne, D.G. 1982, MNRAS, 199, 883 

Blinnikov, S. I., Dunina-Barkovskaya, N. V., & Nadyozhin, D. K. 

1996, ApJ, 106, 171 
Cannizzo, J. K., Lee, H. M. & Goodman, J., 1990, ApJ, 351, 38 
Chen, W. X. & Beloborodov, A. M. 2007, ApJ, 657, 383 
Chincarini, G. et al. 2007, ApJ, 671, 1903 
Dai, Z.G., Wang, X.Y., Wu, X.F. 2006, Sci., 311, 1127 
Di Matteo, T., Perna, T., & Narayan, R. 2002, ApJ, 579, 706 
Fan, Y.Z., Zhang, B. & Proga, D. 2005, ApJ, 635, L129 
Fox, D. B. & Meszaros, P. 2006, New J. Phys., 8, 199 
Gehrels, N., Cannizzo, J.K., & Norris, J. P. 2007, New J. Phys., 9, 

37 

Gehrels, N., et al. 2006, Nature, 444, 1044 

Gorosabel, et al. 2005, A&A, 444, 711 

Itoh, N., et al. 1996, ApJS, 102, 411 Itoh et al. (1996) 

King, A. R., O'Brien, P. T., Goad, M. R., Osborne, J., Olsson, E., 

Page, K. 2005, ApJ, 630, L113 
Kluzniak, W., Ruderman, M. 1998, ApJ, 505, L113 
Kohri, K., & Mineshige, S. 2002, ApJ, 577, 311 
Lattimer, J. M., & Schramm, D. N. 1974, ApJ, 192, L145 
Langanke, K., & Martmez-Pinedo, G. 2001, At. Data Nucl Data 

Tables, 79, 1 

Lazzati, D., Perna, R. & Begelman, M. C. 2008, MNRAS, 388, L15 

Lee, W. H. 2001, MNRAS, 328, 583 

Lee, W.H. & Ramirez-Ruiz, E. 2002, ApJ, 577, 893 

Lee, W. H., Ramirez-Ruiz, E. & Granot, J. 2005, ApJ, 630, L165 

Lee, W. H., & Ramirez-Ruiz, E. 2007, New J. Phys., 9, 17 

Lee, W. H., Ramirez-Ruiz, E. & Page, D. 2005, ApJ, 632, 421 



and short GRBs (see e.g.. IChincarini et al.l ()2007D ). and 
the current proposal applies only to their production in 
short eve nts, given the assu mptions about the progenitor. 
However iKing et al.l ()2005| ) pointed out that fragmenta- 
tion of a rapidly rotating core in the coUapsar model 
could in principle produce secondary episodes of activity 
after the prompt phase. Thus, there could be a common 
underlying cause in both classes in terms of renewed ac- 
tivity in the central engine through large scale flow dy- 
namics. 



We thank J. Cannizzo, N. Gehrels, D. Lazzati, B. Met- 
zger, A. Piro, E. Quataert, M. Rees, L. Roberts, E. 
Rossi and S. Rosswog for comments and discussions, as 
well as the referee for comments and suggestions on the 
manuscript. WL is particularly grateful to A.R. King 
for his glimpses of the obvious. This work was sup- 
ported in part by CONACyT (45845E, WL), PAPIIT- 
UNAM (IN113007, WL), the Packard Foundation (ER- 
R), NASA (Swift NX07AE98G, ER-R) and DOE Sci- 
DAC (DE-FC02-01ER41176, ER-R). DL-C is supported 
by a CONACyT graduate scholarship. WL thanks the 
Department of Astronomy and Astrophysics at the Uni- 
versity of California, Santa Cruz for hospitality. 



Lopez-Camara, D., Lee, W. H. & Ramirez-Ruiz, E. 2009, ApJ, 692, 
804 

MacFadyen, A.L, Woosley, S.E. 1999, ApJ, 524, 262 

Metzger, B. D., Piro, A. L., Quataert, E. 2008, MNRAS, 390, 781 

Monaghan, J.J. 1992, ARA&A, 30, 543 

Nakar, E., 2007, Phys. Rep., 442, 166 

Narayan, R., Piran, T., & Kumar, P. 2001, ApJ, 557, 949 
Nousek, J. A. et al. 2006, ApJ, 642, 389 

O'Brien, P.T., WiUingale, R., Osborne, J. P. 2007, New J. Phys., 8, 
121 

Pacznyski, B., & Wiita, P. J. 1980, A&A, 88, 23 

Perna, R., Armitage, P.J., Zhang, B. 2006, MNRAS, 636, L29 

Popham, R., Woosley, S. E., & Fryer, C. 1999, ApJ, 518, 356 

Proga, D. & Zhang, B. 2006, MNRAS, 370, L61 

Ramirez-Ruiz, E., Merloni, A., Rees, M.J. 2001, MNRAS, 324, 1147 

Rasio, F. A. & Shapiro, S. L. 1994, ApJ, 432, 242 

Rossi, E., & Begelman, M.C., 2009, MNRAS, 392, 1451 

Rosswog, S. 2005, ApJ, 634, 1202 

Rosswog, S. 2007, MNRAS, 376, L48 

Rosswog. S. 2007, RMxAC, 27, 57 

Rosswog, S. & Ramirez-Ruiz, E. 2002, MNRAS, 336, L7 
Setiawan, S., RufTert, M., Janka, H.-Th. 2004, MNRAS, 352, 753 
Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 
Woosley S.E., 1993, ApJ, 405, 273 
Woosley S.E., & Bloom, J.S., 2006, ARA&A, 44, 507 

Zhang, B. et al. 2006, ApJ, 642, 354 

Zhang, B. et al. 2009, Apjsubmitted. larXiv:0902.2419l /l 



